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Abstract 

We consider a simple model of a closed economic system where the total money is conserved 
and the number of economic agents is fixed. Analogous to statistical systems in equilibrium, 
money and the average money per economic agent are equivalent to energy and tempera- 
ture, respectively. We investigate the effect of the saving propensity of the agents on the 
stationary or equilibrium probability distribution of money. When the agents do not save, 
the equilibrium money distribution becomes the usual Gibb's distribution, characteristic of 
non-interacting agents. However with saving, even for individual self-interest, the dynam- 
ics becomes cooperative and the resulting asymmetric Gaussian-like stationary distribution 
acquires global ordering properties. Intriguing singularities are observed in the stationary 
money distribution in the market, as functions of the marginal saving propensity of the 
agents. 
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1 Introduction 

The interacting dynamical nature of any economic sector composed of many cooperatively 
interacting agents, has many features in common with the statistical physics of interact- 
ing systems. In fact, economists like Pareto investigated the power law properties of the 
wealth distributions more than a century ago and Stigler studied the market fluctuations 
by employing Monte Carlo methods more than thirty years back [1]. Motivated by these 
investigations, many efforts are being made recently to apply the statistical physics methods 
to various economic problems. A major part of the recent efforts has gone in investigat- 
ing the nature of fluctuations and their distributions in the stock markets [2]. Wc believe, 
however, a fundamental dynamics occurs in the money market, which affects strongly the 
dynamics of other sectors in the economy. An understanding of the statistical mechanics 
of the money market is therefore essential, and some of its features are very intriguing [3, 
4]. Dragulescu and Yakovenko [5] have shown very recently that for any arbitrary and ran- 
dom sharing but locally conserving money transaction between any two agents in a market, 
the money distribution goes to the equilibrium Gibb's distribution of statistical mechanics: 
P{m) — (l/T) e:K.p{—m/T) where T — M/N, the average money per agent in the market (M 
is the total money of N agents in the market) . This equilibrium distribution is extremely ro- 
bust and various kinds of inter-agent monetary transactions, which locally conserve the total 
money, lead to the Gibb's probability distribution satisfying P{mi)P{m2) = P{mi + 1112). A 
major achievement of this study [5] has been the precise identiflcation of the temperature T 
as the average money per agent in the market. This is due to the fact that the probability 
distribution is normalized and the total money is conserved. This precise identification had 
been missing in many of the previous attempts [3, 4] though the identification from fiscal 
policy considerations was indeed very close [3]. They had considered a market consisting of 
N agents, where initially each one gets any arbitrary share of the total money M in the 
market {J2iLinT'i — M; M and N fixed). The "trade" dynamics goes as follows. Select any 
two arbitrary agents i and j with money and rrij, respectively. These two agents then 
exchange their money through some trade, keeping their total amount of money rrii + rrij 
conserved and no debt is allowed {rrii > at any stage of the trade). There is no other 
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restriction in the trade. Extensive numerical simulations show that this and various modifi- 
cations of trade, like multi-agent transactions, etc., all lead to the robust Gibb's distribution, 
independent of the initial distribution the market starts with. So, most of the agents end-up 
in this market with very little money! Supply of more money in the market (increasing T) 
can increase the width of the distribution, but the most probable money for any agent in 
the market remains zero. It may be mentioned that Ispolatov et al. had studied earlier the 
(non-equilibrium) wealth distributions in various asset exchange models [6] where the trade 
dynamics do not have time-reversal symmetry. Of course, the study of Levy and Solomon 
[7] indicates a power law distribution for the wealth in realistic (possibly non-equilibrium) 
markets or societies. 

Although the model of Dragulescu and Yakovenko [5] is a simple and interesting one to 
start with, a very natural ingredient for any realistic economic agent is missing in the model: 
no economic agent trades with the entire money he or she possesses without saving a part 
of it; saving propensity is too natural a tendency for any economic agent [8]. This saving 
propensity of course varies from agent to agent and even with the accumulation of wealth of 
a single agent. There are also country- wise variations of this saving propensity. We note that 
the fraction of savings A, called the "marginal propensity to save" by economists [8], remains 
fairly constant, independent of the agents. We have taken it to be a constant in the model 
considered here. We show that in presence of this saving factor 0) the money market 
becomes truly cooperative in nature and "critical behaviour" [9] sets in. The multiplicative 
property of the Gibb's distribution P{m), discussed above, was responsible for the absence 
of any cooperative feature of the statistics for A = 0. Once the local or individual measure of 
saving propensity is introduced (A ^ 0), a global order emerges in the entire money market, 
giving a non-vanishing most probable money for each agent in the market and nontrivial 
critical behaviour of the resulting statistics in the equilibrium. 
2 Model and simulation results for money distribution 

We again consider a simple model of the closed economic system where the total amount 
of money M is conserved and the number of economic agents N is fixed. Each economic 
agent i, which may be an individual or a corporate, possesses money mj. An economic agent 
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can exchange money with any other agent through some trade, keeping the total amount of 
money of both the agents conserved. We assume that each economic agent saves a fraction 
A of its money rrii before trading. We again assume that an agent's money must always be 
non-negative and therefore no debt is permitted. Let us now consider that an arbitrary pair 
of agents i and j get engaged in a trade so that their money rrii and rrij change by amounts 
Ami and Arrij to become and m'j, where Ami is a random fraction of (1 — X){mi + mj) 
and Amj is the rest of it. Conservation of the total money in each trade is ensured, as 
earlier. 

We performed computer simulations with fixed number of agents and total money 
M — NT. Most of our simulation results are for — 500 and T — 100. However, we 
checked results for different N values (250 < N < 1000) to check finite size effects, etc., and 
with different T values (upto 10000). Initially we divided the total money M amongst N 
agents equally so that mi = M/N = T for all i. We chose a fixed value of A between zero 
and unity and investigated its effects on the equilibrium distribution of money P{m) in the 
market, giving the (normalized) number of agents P with money m. We choose randomly 
two agents i and j having money and m^, respectively. Then Ami — ^{^ ~ X){'>Tii + mj) 
and Amj = (1 — e)(l — X){mi + mj), where e is a random number between zero and unity. 
Then m' = Am, + Ami and m'j = Xmj + Amj after the trade. Alternatively, this trade 
can also be viewed as m^ m/, mj mj' where m/ = m^ — Am, mj' = m.j + Am 
with Am = (1 — \)[mi — e{mi + mj)]. This can be checked by straightforward substitution. 
These trades were repeated for large number of arbitrary choices of pairs of agents; each of 
these random choices of pairs is considered as one trade and consequently one time unit. The 
typical time upto which we run this algorithm is above 5000 (sometimes upto 50000 for larger 
A). The probability distribution was determined after every 50 time steps till a stationary 
distribution was obtained. We then took an average over 2000 such stationary distributions 
to obtain a smooth distribution. We checked that the equilibrium distribution obtained is 
again extremely robust and does not depend at all on the initial distribution of money in the 
market. We investigated the nature of the equilibrium distribution P(m) for various values of 
A (0 < A < 1). Apart from the stationary distribution, we also investigated the "relaxation" 
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behaviour [9] of the distributions and obtained the time variations of Pi = P{m — T) till 
a steady behaviour was found. We define the relaxation time tr as the earliest time where 
Pi(t) becomes practically independent of time. 

The results for the equilibrium distribution P[m) are shown in Fig. 1, for some values of 
A. The inset shows that the equilibrium distribution is independent of the market size N and 
the average money in the market or temperature T. The real money exchanged randomly in 
any trade is less than the total money, because of the saving by each agent. This destroys the 
multiplicative property of the distribution P(m) (seen earlier for A = 0) and P{m) changes 
from the Gibb's form to the asymmetric Gaussian-like form as soon as a finite A is introduced. 
The A = case was practically a random noise dominating one and therefore effectively a 
non- interacting market. Introduction of a finite amount of saving (A ^ 0), dictated by 
individual self-interest, immediately makes the money dynamics cooperative and the global 
ordering (in the distribution) is achieved. This kind of self-organization in the money market, 
coming out of pure self-interest of each agent, is reminiscent of the "invisible hand" effect 
[8, 10] in the "free market" suggested originally by Adam Smith in 1776. 

The relaxation behaviour of the distribution is shown in Fig. 2, where the time variation 
of Pi{t) for different values of A is shown. Since rrii/T — 1 for all i to start with, Pi{t) starts 
falling from unity in all the cases (not shown in Fig. 2). After some initial rapid decay, we 
see an extremely slow spin-glass type [4, 9] (In t) relaxation behaviour. The inset shows the 
typical variation of the relaxation time tr with A. The dynamics obviously becomes slower 
with increasing A and tr seems to diverge as A approaches unity. Precisely at A = 1, the 
dynamics of course stops. 

An important feature of this humped distribution P{m) at any non- vanishing A is the 
variation of the most probable money mp{\) ( where P{m) becomes maximum) of the agents. 
As discussed before, rup = for A = (Gibb's distribution) and most of the economic agents 
in the market end-up losing most of their money. However, even with pure self-interest of 
each agent for saving a factor A of its own money at any trade, a global feature emerges: 
the entire market ends-up with a most-probable money mp(A). This mp(A) shifts in an 
interesting manner from rUp = (for A = 0) to rUp — > T (for A — > 1). We find that initially. 
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for small A, mp(A) varies very closely as X^/'^ and then it crosses over at A = Ac i:^; 0.45 to A^/^ 
variation. This is shown in Fig. 3, where the curves for A^/^ and A^/^ are also indicated. We 
do not have any idea why this crossover in the exponent for A occurs at any finite A (~ 0.45). 
Also the exponents 1/2 and 1/3 arc very intriguing. Wc checked that this crossover point is 
perhaps also the point where the probability of an agent having zero money Pq = P{m — 
just disappears, as shown in the inset of Fig. 3. The half-width Amp and the peak height 
P^^ of the equihbrium distribution, which scale practically as (1 — A)^/^ and (1 — A)~^/^ 
respectively (see Fig. 4), do not have any irregularity there. In fact, no other property of the 
equilibrium distribution P[m) has any irregularity at Ac. It may be mentioned that Ispolatov 
et al. also observed a singularity in the power law of non-equilibrium growth distribution of 
wealth in their multiplicative asset exchange model having broken time-reversal symmetric 
trade dynamics. 

We noted that although the total distribution assumes some global cooperative feature, 
each individual's money rrii fluctuates randomly. Fig. 5 shows the time variation of the 
money of an arbitrarily chosen individual in the market for two different values of A. The 
inset shows the variation of the time-averaged money (mj) of the agent and its fluctuation 
Amj = {{rrii — {rrii))^) with A after relaxation {t > tr). Since the total money is conserved, 
(rrii) remains constant (= T) here, while Am, goes down with A as (1 — A). This is because 
at any time the agents keep a fixed fraction of their individual money and receive a random 
fraction of the money traded proportional to (1 — A). 
3 Summary and discussion 

We thus considered here a very simple model of money market where the total money M and 
the number of agents is fixed. Each pair of arbitrarily chosen agents in the market exchange 
money with each other through a trade (each trade considered as one time unit). During the 
trade, each agent saves a fraction A of its own money at that time and exchanges randomly 
out of the remaining money, conserving the total amount of money and not allowing any 
debt. We find that for A = 0, the market effectively becomes non-interacting and, no matter 
what initial distribution the agents start with, the resulting money distribution becomes 
the equilibrium Gibb's distribution [5], where most of the agents in the market end-up 
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with very little money. For any non- vanishing A, the equilibrium distribution becomes the 
asymmetric Gaussian-like with the most probable money rup (corresponding to the peak 
in P(jii)) shifting away from m = with increasing A. This global feature, coming out 
of the individual self-interest of saving a part of its own money, can be considered to be 
a demonstration of the self-organization in the market suggested long ago by Adam Smith 
("invisible hand" [8, 10]). Apart from this we find intriguing singularities appearing in the 
equilibrium distribution P{m): nip ~ A^ where v — 1/2 for A < Ac (~ 0.45) and v — 1/3 for 
Ac < A < 1. Also, Amp ~ (1 - Xfl^ and P^^ ~ (1 - A)-V2. 

These singularities in the equilibrium distribution come obviously from the cooperative 
nature of the market interactions induced by the saving propensity of the agents. It may 
be mentioned that while such singular behaviour in the equilibrium money distribution is 
very natural here, in the corresponding physical (gas) system of Newtonian particles one 
gets regular distributions, e.g., the Gibb's distribution (or for that matter, Bose or Fermi 
distributions for quantum particles) and never any singularity, because of the absence of any 
physical equivalent of the "saving" factor there. We thus believe, while regular distributions 
are common for minimally interacting physical many-body systems in equilibrium, singular 
distributions are typical of any working model of the markets. 
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Figure captions 



Fig. 1 : The stationary money distribution P{m) versus money m for different saving 
propensity factor X (N — 500, T — 100). The inset shows the same for different N and T at 
two fixed values of A. The superposition of data indicates the absence of finite size effects 
and the temperature independence. 

Fig. 2 : Relaxation of the distribution: Pi{t) versus the time t (number of trades) for 
different A. Because of extremely slow relaxation, we show in the In t scale for the later 
stage only. The vertical arrows indicate the relaxation time tr. The inset shows the typical 
variation of tr with A; tr diverges as A — > 1. 

Fig. 3 : The variation of the most probable money rrip of an agent (peak position of the 
equilibrium distribution P{m)) as a function of A. The dotted and dashed curves correspond 
to A^/^ and A^/^ respectively. The crossover point (Ac) is indicated by the vertical arrow. 
The inset shows that the probability Pq of agents with zero money also first disappears at 
almost the same point (Ac) indicated by the vertical arrow. 

Fig. 4 : The variation of the distribution half-width Amp as function of A. The dashed 



curve corresponds to a/2. 0(1 — A). The inset shows the peak height Pmj, as function of A. 



The dotted curve here corresponds to 1/^2.4(1 — A). 
Fig. 5 : Time variation of money of an individual agent for two different values of A. 
The inset shows the variations of (mj) and its fluctuations Am^ as functions of A. 
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